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Using BABMi measurements of the inclusive electron spectrum in B ^ X^ev decays and the inclusive 
photon spectrum in B ^ Xs'y decays, we extract the magnitude of the CKM matrix element Vub- 
The extraction is based on theoretical calculations designed to reduce the theoretical uncertainties 
by exploiting the assumption that the leading shape functions are the same for all 6 ^ g transitions 
(g is a light quark) . The results agree well with the previous analysis, have indeed smaller theoretical 
errors, but are presently limited by the knowledge of the photon spectrum and the experimental 
errors on the lepton spectrum. 

PACS numbers: 13.20.He, 12.15.Hh, 12.38.Qk, 14.40.Nd 



I. INTRODUCTION 

The determination of the magnitude of the Cabibbo- 
Kobayashi-Maskawa matrix element Vub from charmless 
semileptonic B decays is complicated by the fact that 
over most of the phase space B Xctv decays domi- 
nate and are very difficult to distinguish from the signal 
B Xuiv decays. Here Xc and Xu refer to hadrons, 
mostly mesons, with and without charm. To reduce im- 
pact of the dominant B XJ.^ background, partial 
rates for B — > Xu^v decays are measured in regions of 
phase space where these background decays are forbidden 
or highly suppressed, for instance, near the endpoint of 
the energy spectrum of the charged lepton I. The partial 
decay rates are extrapolated to the total rate by compar- 
ing the experimentally measured rate with a theoretical 
prediction. 

Theoretically, the most precise predictions can be 
made for the total B Xu^,v decay rate. Accounting 
for restriction in phase space is difhcult because decay 
spectra close to the kinematic limit are susceptible to 
non-perturbative strong- interaction effects. Theoretical 
tools for the calculations of the partial inclusive decay 
rates are QCD factorization and local operator prod- 
uct expansions (OPE). These calculations separate non- 
perturbative from perturbative quantities, and use ex- 
pansions in inverse powers of the 6-quark mass, mf,, and 
in powers of the strong coupling as [1] . At leading or- 
der in Aqcd/w{,, the non-perturbative bound-state ef- 
fects are accounted for by a shape function describing 
the "Fermi-motion" of the b quark inside the B me- 
son. These shape functions cannot be calculated. Dif- 
ferent shapes have been proposed, and parameters defin- 
ing these shapes have to be extracted from data. This 



introduces significant additional hadronic uncertainties. 
At leading order, these shape functions are assumed to 
be universal functions for & — > g transitions, where g is a 
light quark, either s or u. 

In the past, we have extracted the non-perturbative pa- 
rameters of these shape functions, the 6-quark mass mf, 
and its mean kinetic energy squared, J (/i) , from the mo- 
ments of the inclusive photon spectrum in B ^ Xg^ de- 
cays, as well as from moments of the hadron mass and the 
lepton energy distributions in _B ^ Xcf-v decay. These 
parameters depend on the choice of the renormalization 
scale, /i. 

It was suggested many years ago [2] that \Vub\ can be 
extracted with smaller theoretical uncertainties by com- 
bining integrals over the lepton spectrum in B — > X^tv 
decays with weighted integrals over the photon spectrum 
from B — > decays, each above a cut-off energy Eq. 
The underlying assumption is that the QCD interactions 
affecting these two processes are the same and thus will 
cancel to first order in the appropriate ratio of weighted 
decay rates. The advantage of this approach is that it re- 
duces the sensitivity to the choice of the shape-function 
parameterization and thus avoids uncertainties that are 
difficult to quantify. 

In the following, we extract \Vub\ using two differ- 
ent prescriptions, one proposed by Leibovich, Low, and 
Rothstein (LLR) [3-5] , and the other based on more re- 
cent calculations by Lange, Neubert, and Paz (LNP) [6, 
7] . The two prescriptions use different calculations of the 
weight functions, and thus result in different estimates of 
the theoretical uncertainties. 

The BaBAR Collaboration was the first to apply the 
LLR prescription [3, 4] to extract \Vub\, based on mea- 
surements of the hadron mass spectrum in charmless 
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FIG. 1: The difFerential branching fraction for B — > Xuev 
decays [9] as a function of the charged lepton momentum, 
measured in the T{4S) rest frame. The error bars indicate 
the statistical errors; the histogram represents the theoretical 
prediction [11]. 



semileptonic B decays and the inclusive photon spec- 
trum [8]. 



rived from the predicted electron spectrum [11], calcu- 
lated with the shape-function parameters determined by 
a global fit [12] to moments of inclusive distributions in 
semileptonic and radiative i?-meson decays. The system- 
atic uncertainties for this transformation are estimated 
by varying the shape-function parameters within their 
experimental uncertainties. The resulting partial branch- 
ing fractions for different values of the electron-energy 
cut-off, Eq, measured in the T(4S') and the B-meson rest 
frames, and correction factors relating the two, are listed 
in Table I. 



TABLE I: Partial branching fractions for B X^ev decays 
in units 10~^, integrated over the energy range from Eq to 
2.6 GeV/c, both in the T(45) rest frame and the B-meson rest 
frame, as well as the correction factor relating the two. The 
uncertainty of the correction factor reflects the uncertainty in 
the assumed shape of the spectrum. 



Eo AI3{Eo) ■ 10^ Correction AB{Eo) ■ 10^ 

[GeV] T(4S) rest frame factor B rest frame 

""2l) 0.572 ± 0.041 ± 0.065 1.002 ± 0.005 0.573 ± 0.077 

2.1 0.392 ± 0.023 ± 0.038 0.994 ± 0.008 0.390 ± 0.044 

2.2 0.243 ±0.011 ±0.020 0.973 ± 0.016 0.236 ± 0.023 

2.3 0.148 ±0.006 ±0.010 0M5t.o '^tl 0.135 ±0.012 

2.4 0.075 ± 0.004 ± 0.006 0.77211^:^81 0.058 ± 0.008 



II. EXPERIMENTAL INPUT 



The experimental inputs for this analysis are the pub- 
lished BaBAR measurements of the inclusive electron en- 
ergy spectrum in _B — > XuCv decays [9], and of the in- 
clusive photon-energy spectrum in B Xgj decays [10]. 
These measurements are based on a data sample cor- 
responding to a total integrated luminosity of about 
80 fb^^. The measured spectra will be integrated above 
an energy Eq, measured in the _B-meson rest frame. 



A. Inclusive Lepton Spectrum in S — > X^ev 
Decays 



B. Inclusive Photon Spectrum in B — >■ Xsf Decays 

The inclusive photon-energy spectrum [10] in B —^ 
Xs'j decays, above 1.9 GeV, measured in the B-meson 
rest frame, is shown in Fig. 2. This spectrum is mea- 
sured as the sum of the photon spectra in 38 exclusive 
decay modes. The data are fully corrected for detection 
efficiencies and energy resolution. They are normalized 
to the total number of charged and neutral _B-meson de- 
cays and presented as differential branching fractions. 



The inclusive electron-energy spectrum above 2.0 GeV, 
measured in the T(4S') rest frame, is shown in Fig.l. The 
data are fully corrected for detection efhciencies as well 
as final state radiation and bremsstrahlung. They are 
normalized to the total number of charged and neutral B- 
meson decays and are presented as differential branching 
fraction. 

For the extraction of \Vub\ we need to transform the 
measured partial branching fraction from the T(4S') rest 
frame to the branching fraction in the B-meson rest 
frame, integrated over the spectrum above an energy 
cut-off at Eq . This is done using correction factors de- 



III. THEORETICAL FRAMEWORK FOR THE 
EXTRACTION OF \Vub\ 

A. Method I by Leibovich, Low, and Rothstein 

A. K. Leibovich, L Low, and 1. Z. Rothstein [3, 4] 
proposed a method for the extraction of the ratio 
iKibl^/jVtbVj* [^ without invoking knowledge of the shape 
function. Their calculation relates | Kibj^/jVtbVj* [^ to the 
experimentally measured differential branching fractions 
for B Xutv and B — > Xg^ decays. 
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FIG. 2; The differential branching fraction of inclusive B 
Xs"/ decays as a function of the photon energy in the B rest 
frame. The error bars indicate the statistical errors. The data 
are compared to the theoretical prediction [11] using shape 
function parameters mt — 4.67 GeV, /ij — 0.16 GeV^ that 
have been obtained from the fit to E-y spectrum [10]. 
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The integration variables are 

xb^2Ei/Mb, ub = 2E^/Mb, 

with limits Xg and x'^/p; Mb refers to the i3- meson 
mass. The weight function is approximately linear as 
a function of ub. 
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FIG. 3: The weight functions for Method I (LLR) for two 
values of = 2.0 GeV and Eq = 2.3 GeV, corresponding to 
different values of p = 0.9983 (solid lines), p = 0.995 (dashed 
lines), and non-resummed weight function (dotted lines). 



The Wilson coefficients cf'\mh) , C^^nib), etc., com- 
puted in NLO, can be found in [13, 14]. The argu- 
ment oi K, y ^ 2;^ln(l - asPolxB/uB)i diverges for 
(^sPqIxb/ub = 1- To avoid this pole, the integration lim- 
its are set to < p ub with p < 0.999. For smaller 
values of p, the weight function deviates and thus the 
extracted value of |yi,f,| changes, and the uncertainty of 
this scheme increases. 

To check the impact of the resummation of the Su- 
dakov logarithms we also tried a non-resummed version 
of the LLR weight function [5] : 



^ dxB(^l-Hl-XB)^ 



^(^-^-^ln(l-^))). (3) 
TT 2 9 9 Ub 



This calculation (Eq. 1 — 3) includes NLO perturbative 
corrections, but does not take into account power cor- 
rections. Estimates of the theoretical uncertainties are 
discussed in Section 4. The weight function is shown in 
Fig. 3, for Eo = 2.0 GeV and -Bo = 2.3 GeV and different 
values of the parameter p. 

M. Neubert [15] performed calculations that are quite 
similar and give results that are very close to the ones 
obtained for this method. 



B. Method II by Lange, Neubert, and Paz 

The second method for the extraction of \Vub\ is an 
application of the more recent two-loop calculations by 
B. Lange, M. Neubert, and G. Paz [6], and by B. 
Lange [7] , performed for the high end of the lepton and 
photon energy spectra in i? — > X^^i/ and B — > Xsj de- 
cays. 
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FIG. 4: The weight function (in units of ps ^) for Method 
II (LNP), for Eo = 2.0 GeV (upper curve) and Eo = 2.3 GeV 
(lower curve). 



Unlike the previously described methods, this one re- 
lates \Vub\ to the measured partial B branch- 
ing fraction and normalized photon spectrum in the 
S — > decay, 



Mb/2 

^ J dEi{dB{B^XJiy)/dEi) 

Eq 

Mb/2 

/ dE^ w{E^, Eo)S{E^) + r,hc(£;o) 

Eq 

1 dr{B ^ X,7) 



T{B ^ X«7)b,>e„,„ 



dE^ 



(4) 



(5) 



where rrhc(£'o) represents residual hadronic power cor- 
rections, which in [6] were absorbed into the weight func- 
tion, ii^niin = 1-90 GeV is chosen as the lower limit for the 
normalization of the photon energy spectrum; the corre- 
sponding branching fraction is 3.95 x lO"'^. 

This calculation contains pcrturbative corrections at 
NNLO at the so-called "jet scale", /Xi ~ a/w^Aqcd, 
and at NLO at the so-called "hard scale", fih ~ Wf,. 
Also included are the first-order power corrections, which 
are separated into two parts: kinematic corrections and 
residual hadronic corrections. The kinematic corrections 
do not introduce hadronic uncertainties and are applied 
directly to the weight function. The residual hadronic 
corrections include subleading shape functions for which 
the functional form is unknown. This introduces signifi- 
cant theoretical uncertainties. The weight functions, cal- 
culated for Eq = 2.0 GeV and Eq = 2.3 GeV, are shown 
in Fig 4. 



IV. DETERMINATION OF \Vub\ 

A. Calculation of weighted integrals of B 
photon spectrum 



For each of the two methods we extract \Vub\ from a 
ratio of the B X^lv partial branching fraction and 
a weighted integral over the photon spectrum for B 
Xgj. Each method introduces a specific weight function 

w{E^i,Eq). The value of the weighted integral of the 
photon spectrum above a minimum energy Eq, 



Mb/2 



dE^w{E^,Eo)dT{E^)/dE^, (6) 



is taken as a sum over bins in the photon energy in the 
B-rest frame, 

i[E,) = Y,y^{E,i,E,)(^^^)AE,,. (7) 

i 

The uncertainty for this sum is estimated using the 
standard error propagation for a linear combination of 

random variables: 

(7^ =Y^w{E^i,Eo)w{E^j,Eo)ilE^AE^jVij, (8) 



where Vij is a covariance matrix of the differential rate 
dr{Ej)/dE^ for different photon energy bins. Assuming 
uncorrelated statistical errors and correlated systematic 
errors with known correlation matrix Rij , the covariance 
matrix is of the form 



''stat'^'stat^ii 



(9) 



The correlation matrix Rij was provided by the BaBAR 
Collaboration [10]. 



B. Results on |Vub| and Error Estimation 

1. Method I (LLR) 

The results of the extraction of IKisl/lVjbVj* | based 
on Method I arc presented in Tabic II. The theoreti- 
cal uncertainties of the NLO calculations have been es- 
timated [3, 4] to be ©(Aqcd/Mb), resulting in a rela- 
tive error of about 6% for |V^&|/|Vt6V^(* |. The extraction 
of iKtfcl/IVtftl^* I was done with rcsummcd Sudakov log- 
arithms. As a cross check, the extraction was also per- 
formed without resummation, resulting in decreases of 
the weight function by ~ 6% at Eq = 2.0 GeV and by up 
to ~ 12% at 2.4 GeV. 

To translate these results into \Vub\ we exploit the con- 
straints of the unitarity of the CKM matrix resulting in 
\Vtb\ = 1 + C(A4), where A = 0.226 is the sine of the 
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TABLE II: The results for |V^u6|/|Vti,V't3| and \Vub\ for Method I (LLR). The extraction uses the weight functions based on 
resummed Sudakov logarithms with p = 0.9983. The first error represents the error from the measured B —^ Xy,ev partial 
branching fraction, the second error is from the measured B Xsj spectrum, and the third is the estimated theoretical 
uncertainty. The fourth error on \Vub\ is from the \Vts\ uncertainty. 



Eq [GeV] \Vub\/\VthVt*s\ IKfl-lO^ 

2.0 0.106 ± 0.007 ±0.007 ±0.006 4.28 ± 0.29 ± 0.29 ± 0.26 ± 0.28 

2.1 0.100 ±0.006 ±0.006 ±0.006 4.06 ± 0.23 ± 0.25 ± 0.24 ± 0.27 

2.2 0.093 ± 0.005 ± 0.005 ± 0.006 3.78 ± 0.18 ± 0.21 ± 0.23 ± 0.25 

2.3 0.091 ± 0.004 ± 0.005 ± 0.006 3.69 ± 0.16 ± 0.19 ± 0.22 ± 0.25 

2.4 0.090 ± 0.006 ± 0.004 ± 0.006 3.64 ± 0.25 ± 0.17 ± 0.22 ± 0.24 



Cabibbo angle. Using \Vts\ = (40.6 ± 2.7) • IQ-^ [16], we 
calculate \Vub\ (see Table II). 

The results show high stability with respect to varia- 
tions of the Icpton energy cnt in the range from 2.0 to 
2.4 GeV. The partial charmless semileptonic branching 
fraction decreases in this range from 25% to 2.3% of the 
total B Xulv branching fraction [9]. The observed 
stability is quite surprising because near the endpoint of 
the lepton energy spectrum the theoretical calculations 
are expected to break down, and this should lead to in- 
creasing theoretical uncertainties. 

2. Method 11 (LNP) 

The results on the extraction of \Vub\ based on Method 
II are presented in Table III. The estimated individual 
theoretical uncertainties [7] are added in quadrature to 
obtain the total theoretical error. The first order non- 
perturbative hadronic power corrections, not considered 
in Method I, are split into two contributions: kinematic 
corrections, which are included in the weight function and 
depend on the scale for the calculation of kinematic power 
corrections, /I, but are independent of intermediate scale, 
fii, the hard scale, ^/j, and residual hadronic corrections. 
These corrections depend on unknown subleading shape 
functions, which affect the B X^lu and B — » Xgj 
spectra in different ways and thus do not cancel in the 
weighted ratio in the integrals. This uncertainty, crhad, is 
estimated as suggested by Lange [7], namely by varying 
the shape and the parameters of the subleading shape 
functions. 

The error Cpert is the uncertainty of the NNLO ap- 
proximation of the shape function and the LO approxi- 
mation of the power corrections. It was estimated from 
the dependence of the sum of the weighted integral over 
the photon spectrum and residual hadronic correction on 
the hard scale, fXh = mb/V^, on the intermediate scale, 
Hi = 1.5 GeV, and on the scale for the power corrections, 
/I = 1.5 GeV. All three scales are varied by factors of \/2 
and l/-\/2 relative to their default values, and the largest 
variation of the weighted integral is taken as the estimate 
of the perturbative uncertainty. 

The errors CTmu and cTnarm are due to uncertainties of 



the parameters that are inputs to the calculation and 
were varied within their stated errors: TOf, = 4.61 ± 0.06 
GeV and mc/rrib = 0.222 ± 0.027, = 90 ± 25 MeV. 
Here irib and /i^ are defined in the shape-function scheme 
at a scale /u* = 1.5 GeV, nig and the ratio rric/mb are 
evaluated in the MS scheme at 1.5 GeV, where the ratio 
is scale-invariant. The imcertainties on A2 = 0.12 GeV^ 
and /x^ = 0.25 ±0.10 GeV^ only enter into the subleading 
shape functions, for which the uncertainties are assessed 
separately. 

The error (Jnorm represents the uncertainty of the nor- 
malization of the photon spectrum, which is estimated to 
be about 6% [7]. 

We observe in Table III a significant increase in the 
extracted value of \Vub\ for Eq > 2.2 GeV, and the 
theoretical uncertainties also show a rapid growth for 
Eo>2.1 GeV, reaching 100% above Eq ~ 2.3 GeV. This 
is due to the power correction ri.hc(-E'o) in the denomina- 
tor of the Eq.4. This correction is negative and almost 
independent of Eq. As £"0 increases, the integral over 
the photon spectrum decreases and contribution from 
the rrhc(i?o) to the total uncertainty increases. This is 
not unexpected, because the effect of non-perturbative 
hadronic corrections increases in the region close to kine- 
matic endpoints for both decays. 



V. CONCLUSION 

We have extracted the CKM matrix element \Vub\ us- 
ing published BaBar measurements of the inclusive lep- 
ton spectrum in B — > X„e^ decays and inclusive photon 
spectrum in _B ^ Xg"/ decays. By using the ratios of 
the weighted spectra for these two decays, the results are 
expected to be less model dependent than previous mea- 
surements relying on the extraction of the shape func- 
tions from data and specific parameterizations of these 
functions. 

For comparison, the results for £^0 = 2.0 GeV for both 
methods are presented in Table IV, together with the 
BaBAR shape-function based measurement [9]. Figure 5 
shows the dependence of the results on the lepton energy 

cut-off, Eq. 

The results agree well within their stated uncertain- 
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TABLE III: The results for \Vub\ based on Method II (LNP). The first error represents the error from the measured B — » Xuev 
partial branching fraction, and the second from the measured B Xs^ spectrum. The contributions to the theoretical error, 
as described in the text, are given together with the total theoretical error. 



Eo [GeV] 


IKftI • 10' 


0"had ■ lO' 


0"pcrt • 10"* 


cr™, • 10=* 


^parm ' 10 


^norm ' 10 


fthcory ' lO' 


2.0 
2.1 
2.2 
2.3 


4.40 ± 0.30 ± 0.41 
4.55 ±0.26 ±0.45 
5.01 ±0.24 ±0.60 
6.99 ±0.31 ±1.60 


+0.08 
-0.07 
+0.15 
-0.13 
+0.42 
-0.33 
+3.90 
-1.44 


+0.03 
-0.01 
+0.11 
-0.05 
+0.40 
-0.21 
+3.02 
-1.00 


+0.13 
-0.12 
+0.16 
-0.15 
+0.25 
-0.22 
+0.75 
-0.58 


+0.01 
-0.00 

±0.01 
±0.01 

+0.04 
-0.03 


±0.17 
±0.21 
±0.32 
±0.92 


+0.23 
-0.22 
+0.32 
-0.29 
+0.71 
-0.55 
+5.07 
-2.05 



TABLE IV: Comparison of the \Vub\ extraction for Eo = 
2.0 GeV for the two methods. The first error refiects the 
uncertainty in the measurements of the B — > Xuli^ lepton 
spectrum, the second error is due to the measurement of the 
B — > Xsy photon spectrum. For the shape-function based 
analysis, the second error originates from the extraction of 
the shape-function parameters, in this case based on both 
the inclusive photon spectrum as well as hadron-mass and 
lepton-energy moments from B Xclv decays. The third is 
the theoretical uncertainty. The fourth error for Method I is 
due to the uncertainty of \Vts\- 



Method 


\Vub\- 10-' 


LLR [3, 4] 
LNP [6, 7] 


4.28 ± 0.29 ± 0.29 ± 0.26 ± 0.28 
4.40 ±0.30 ±0.41 ±0.23 


SF-based analysis [9] 


4.44 ± 0.25 ^ols ± 0.22 




FIG. 5: Comparison of \Vub\ values extracted from the two 
different calculations as a function of the lepton energy cut- 
off, Eq. The errors bars represent the experimental and the 
total error. 



ties. Method's I results in values of \ Vub\ and errors that 
appear to be unaffected by the restriction of the phase 
space near the kinematic limit. The theoretical errors are 
somewhat smaller than for Method II, and also smaller 
than for the previous extraction method, but there is an 



additional error due to the uncertainty of | Vts \ ■ 

The values of \Vub\ extracted by using Method II for 
the lepton energy cut-off at Eq — 2.0 GeV and 2.1 GeV 
agree well with those of Method I. The theoretical error 
contributions are estimated in detail. The smallest total 
theoretical error is obtained for Eq — 2.0 GeV; it is about 
5%. However, at larger Eq^ there are several increasing 
contributions to the theoretical uncertainties, in particu- 
lar the first-order non-perturbative hadronic power cor- 
rections that are not dealt with in Method I. This is why 
the error contribution from the weighted integral over 
the photon spectrum increases very rapidly as the phase 
space gets more restricted. 

At present, for i^o = 2.0 GeV, the experimental errors 
on both the B X^dy branching fraction and the inte- 
gral over B Xg^ are about 12%; this means that there 
are opportunities to improve the accuracy of | Vub \ signif- 
icantly by reducing the experimental errors of these two 
spectra with more data available now and in the future. 
Also, extending i?o to 1.9 GeV or lower, would allow us 
to establish more clearly the stability of the results in 
this region. 

It would be of interest to perform a similar analysis for 
the hadron-mass distribution in i? — s- Xuiiy decays, with 
increasingly tighter restrictions on the mass of and 
assess whether this measurement is indeed more robust 
than the measurement of the lepton spectrum near the 
kinematic endpoint. 
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